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Abstract 

The evolution of the dark energy parameter within the scope of a spatially homogeneous and isotropic 
Friedmann-Robertson- Walker (FRW) model filled with perfect fluid and dark energy components is studied 
by revisiting the recent results (Amirhashchi et al. in Int. J. Theor. Phys. 50: 3529, 2011). The two sources 
are claimed to interact minimally so that their energy momentum tensors are conserved separately. To prevail 
the deterministic solution we consider a time dependent deceleration parameter (DP) i.e. q — — b(a(t)), 

which yields a scale factor a = [sinh(crf)] ~ , where a and n are arbitrary constants. This provides for a model 
which generates a transition of the universe from the early decelerating phase to the recent accelerating 
1 phase. It is observed that the transition red shift (z t ) for our derived model with qo = —0.73 is = 0.32. 

6J[) This is in good agreement with the cosmological observations in the literature. The physical aspects of the 

corresponding solutions are also discussed. 
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1 Introduction 



The 2011 Nobel Prize in Physics has been awarded to Perlmutter, Riess and Schmidt for their very careful and 
independent observations of distant supernovae, apparently showing the rate of expansion of the Universe is 
accelerating. In 1998, the promulgated observations of Type la supernova (SNela) established that our universe 
is presently accelerating (Perlmutter et al. [HE]; Riess et al. [3]) and recent observations of SNela of high 
sureness level (Tonry et al. [4]; Riess et al. [5]; Clocchiatti et al. [6]) have further confirmed this. In addition, 
measurements of the cosmic microwave background (CMB) (Bennett et al. [7. j) and large scale structure (LSS) 
(Tegmark et al. [5]) strongly indicate that our universe is dominated by a component with negative pressure, 
dubbed as dark energy. These also indicate that the universe has a flat geometry on large scales. Because there 
is not enough matter in the universe — ether ordinary or dark matter — to produce this flatness, the difference 
must be attributed to a "dark energy" . This same dark energy causes the acceleration of the expansion of the 
universe. In addition, the effect of dark energy seems to vary, with the expansion of the universe slowing down 
and speeding up over different time. The Wilkinson Microwave Anisotropy Probe (WMAP) satellite experiment 
suggests 73% content of the universe in the form of dark energy, 23% in the form of non-baryonic dark matter 
and the rest 4% in the form of the usual baryonic matter as well as radiation. 



High-precision measurements of expansion of the universe are required to understand how the expansion rate 
changes over time. In general relativity, the evolution of the expansion rate is parameterized by the cosmological 
equation of state (the relationship between temperature, pressure, and combined matter, energy, and vacuum 
energy density for any region of space) . Measuring the equation of state for dark energy is one of the biggest 
efforts in observational cosmology today. The DE model has been characterized in a conventional manner by 
the equation of state (EoS) parameter w de = which is not necessarily constant, where p de is the energy 

density and p de is the fluid pressure (Carroll and Hoffman [5]). The oj de lies close to —1: it would be equal 
to —1 (standard ACDM cosmology), a little bit upper than —1 (the quintessence dark energy) or less than — 1 
(phantom dark energy). While the possibility us <C — 1 is ruled out by current cosmological data from SN la 
(Supernovae Legacy Survey, Gold sample of Hubble Space Telescope) (Riess et al. [5]; Astier et al. [IQ]), CMB 
(WMAP, BOOMERANGE) (Eisentein et al. [Ell; MacTavish et al. [T2]) and large scale structure (Sloan Digital 
Sky Survey) data (Komatsu et al. [13]), the dynamically evolving DE crossing the phantom divide line (PDL) 
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(uj = — 1) is mildly favoured. The simplest candidate for the dark energy is a cosmological constant A, which has 
pressure p( de ^ = —p( de \ Specifically, a reliable model should explain why the present amount of the dark energy 
is so small compared with the fundamental scale (fine-tuning problem) and why it is comparable with the critical 
density today (coincidence problem) (Copeland et al. [II])- That is why, the different forms of dynamically 
changing DE with an effective equation of state (EoS), a/ de ) = p( de )/p( de ) < —1/3, have been proposed in the 
literature. Some other limits obtained from observational results coming from SNe la data (Knop et al. |15j ) 
and combination of SNe la data with CMBR anisotropy and galaxy clustering statistics (Tegmark et al. [5]) 
arc —1.67 < uj < —0.62 and —1.33 < uj < —0.79, respectively. The recent results in 2009, obtained after a 
combination of cosmological datasets coming from CMB anisotropics, luminosity distances of high redshift type 
la supernovae and galaxy clustering, constrain the dark energy EoS to —1.44 < uj < —0.92 at 68% confidence 
level (Komatsu et al. [T3]; Hinshaw et al. [TB"]). 

Recently, Amirhashchi et al. [T7II18] . Pradhan et al. |19j. Saha et al. |20j have studied an interacting and non- 
interacting two-fluid scenario for dark energy models in FRW universe. Kumar [21] studied some FRW models of 
accelerating universe with dark energy. In this paper we study the evolution of the dark energy parameter within 
the framework of a FRW cosmological model filled with two fluids (barotropic and dark energy) by revisiting 
the recent work of Amirhashchi et al. (2011b) and obtained more general results. The cosmological implications 
of this two-fluid scenario will be discussed in detail in this paper. In doing so the two sources are claimed to 
interact minimally so that their energy momentum tensors are conserved separately. The out line of the paper is 
as follows: In Sect. 2, the metric and the basic equations are described. Sections 3 deals with the solutions of 
the field equations. The physical significances of the model s are discussed in Sect. 4. Finally, conclusions are 
summarized in the last Sect. 5. 



2 Model and field equations 

In standard spherical coordinates (x 1 ) — (t, r, 6,cp), a spatially homogeneous and isotropic FRW line-element has 
the form (in units c = 1) 



ds 2 = -dt 2 +a 2 (t) 



dr 2 



1 — kr 2 



+ r 2 (d9 2 + 



2 i -2 

1 sin 



(1) 



where a(t) is the cosmic scale factor, which describes how the distances (scales) change in an expanding or 
contracting universe, and is related to the redshift of the 3-space; k is the curvature parameter, which describes 
the geometry of the spatial section of space-time with closed, flat and open universes corresponding to k = — 1, 
0, 1, respectively. The coordinates r, 9 and <f> in the metric ([1]) are comoving coordinates. The FRW models 
have been remarkably successful in describing the observed nature of universe. 



The Einstein's field equations in case of a mixture of perfect fluid and DE components, in the units 8wG = 
c = 1, read as 

Ugi = Ti, (2) 



where T- = T^" 1 ^ -\-T^ de ^ is the overall energy momentum tensor with T^ m ^ and T^'^ as the energy momentum 
tensors of ordinary matter and DE, respectively. These are given by 



T^ 1 = diag 



-n( m ) r)( m ) r>(" 1 ) 

f 3 y i y ) y 



and 



= diag 

T<f e)i = diag 
= diag 



1, uj {m \ w (m) , w (m) 



.pCdO, p(de) > p (de) i p (de) 



l,uj( de \uj( de \uj^ 



(3) 



(4) 



where p^ m ^ and p( m ) are, respectively the energy density and the isotropic pressure of the perfect fluid compo- 
nent or ordinary baryonic matter while u/ m ) = p( m ) /p^" 1 - 1 is its EoS parameter. Similarly, p( rfe ) and p( de ) are, 
respectively the energy density and pressure of the DE component while 

Jde) = p {de) 

I p( de ) is the corresponding 
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EoS parameter. 

In a comoving coordinate system, the field equations ((2J, for the FRW space-time ([1]), with ([3]) and (UJ, read 



as 

. a a 



2^ + " + * = - w H p (") _ w (*0 p (*0 , ( 5 ) 
a a- 2 

3^ + 4")-P (m) +P (de) - (6) 



Here the over dot denotes derivative with respect to t. 

The law of energy conservation equation T l] .- = yields 

p (m) + 3(1 + u^)p^H + p {de) + 3(1 + u^)p^H = 0, (7) 
where H = - is the Hubble parameter. 

3 Solution of field equations 

The field equations ([5]) and ([BJ involve five unknown variables, viz., a, aj( m \ to^ de \ p^ and p( de \ Therefore, 
to find a deterministic solution of the equations, we need three suitable assumptions connecting the unknown 
variables. 

In order to solve the field equations completely, we first assume that the perfect fluid and DE components 
interact minimally. Therefore, the energy momentum tensors of the two sources may be conserved separately. 

Following Akarsu and Kilinc [22], first we assume that the perfect fluid and DE components interact mini- 
mally. Therefore, the energy momentum tensors of the two sources may be conserved separately. 

The energy conservation equation (T,^ 1 ^ = 0) of the perfect fluid leads to 

p (m) +3(l + cj( m V m )ff = 0, (8) 
whereas the energy conservation equation (T.^. *^ = 0) of the DE component yields 

p (de) +3(1 + J d ^)p {d ^H = 0. (9) 



Following Akarsu and Kilinc [22, 23, 24 , and Kumar and Yadav [25], we assume that the EoS parameter of the 
perfect fluid to be a constant, that is, 

W <™> = ^ = const., (10) 

while u/ de ) has been admitted to be a function of time since the current cosmological data from SN la, CMBR 
and large scale structures mildly favor dynamically evolving DE crossing the PDL as discussed in previous section. 

Eq. © can be integrated to lead 

p^=p a-^ m) \ (11) 
where po is a positive constant of integration. 



Firstly, we define the deceleration parameter q as 



(HI 



H + H 



2 



The motivation to choose such time dependent DP is behind the fact that the universe has accelerated expansion 
at present as observed in recent observations of Type la supernova (Riess et al. 0; Perlmutter et al. pQ; Tonry 
et al. [3]; Riess et al. [5]; Clocchiatti et al. |lj) and CMB anisotropies (Bennett et al. 0; de Bernardis et al. 
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[26] ; Hanany et al. ^ZT\), and decelerated expansion in the past. Also, the transition redshift from decelerated 
expansion to accelerating expansion is about 0.5. Now for a Universe which was decelerating in past and acceler- 
ating at the present time, the DP must show signature flipping (see the Refs. Padmanabhan and Roychowdhury 
[28] ; Amendola [35] ; Riess et al. [3D] )■ So, in general, the DP is not a constant but time variable. 



Equation (|f 21) may be rewritten as 



-+&^ = 0. (13) 
a a^ 



In order to solve the Eq. (|16|) . we assume 6 = 6(a). It is important to note here that one can assume 
6 = b(t) = b(a(t)), as a is also a time dependent function. It can be done only if there is a one to one 
correspondences between t and a. But this is only possible when one avoid singularity like big bang or big rip 
because both t and a are increasing functions. 

The general solution of Eq. (fTB"]) with the assumption 6 = 6(a), is obtained as 

J ef« da da = t + k, (14) 

where k is an integrating constant. 

One cannot solve (|14[) in general as 6 is variable. So, in order to solve the problem completely, we have to 
choose J ^ da in such a manner that (fT4"|) be integrable without any loss of generality. Hence we consider 

J h -da = \nf{a), (15) 

which does not affect the nature of generality of solution. Hence from (fi"4|) and (fT5|) , we obtain 

f{a)da = t + k. (16) 



Of course the choice of /(a), in (|16| . is quite arbitrary but, since we are looking for physically viable models of 
the universe consistent with observations, we consider 

/(a) = " <t "~ 1 , (17) 
V 1 ay/1 + a 2 ™ 

where a is an arbitrary constant and n is a positive constant. In this case, on integrating Eq. (|16[) and neglecting 
the integration constant fc, we obtain the exact solution as 

a(t) = (sinh(at))». (18) 

This relation (fT5)l generalizes the value of scale factor obtained by Amirhashchi et al. [T5] and Pradhan et al. 
[31] in connection with the study of dark energy models respectively in FRW and Bianchi type-Uio space-times. 

From (|18p . we obtain the time varying deceleration parameter as 

= -^ = n [l_(tanh(crf)) 2 l -1. (19) 

From Eq. (|i"9"]). we observe that q > for t < itanh _1 (l - ~)2 and q < for t > itanh _1 (l - It is 

also observed that for < n < 1, our model is in accelerating phase but for n > 1, our model is evolving from 
decelerating phase to accelerating phase. Also, recent observations of SNe la, expose that the present universe 
is accelerating and the value of DP lies to some place in the range — 1 < q < 0. It follows that in our derived 
model, one can choose the value of DP consistent with the observations. Figure 1 depicts the variation of the 
deceleration parameter (g) versus time (t) which gives the behavior of q for different values of n. It is also clear 
from the figure that for n < 1, the model is evolving only in accelerating phase whereas for n > 1 the model is 
evolving from the early decelerated phase to the present accelerating phase. 



Using (fT5|) , the model ([T]) becomes 

2 T dr 2 
ds 2 = -dt 2 + (sinh(crf))^ 



1 — fcr 2 



(d9 2 + sin 2 8dcj) 



(20) 
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Figure 1: The plot of deceleration parameter q vs time t. Here a = 1. 

4 Physical and Geometric properties of the Model 

The spatial volume (V), Hubble parameter (H), expansion scalar (9), energy density (p*" 1 -*) of perfect fluid, DE 
density (p( de " > ) and EoS parameter (oj^ de " > ) of DE, for the model (|2T)1) are found to be 



V = (sinh(ai))", 

9 = 3H = — coth(at), 
n 

(m) _ PO 



(sinh(ai))" (1+ " (m,) ' 
3k 



Po 



Jde) 



(sinh(at))^ (sinh(at))t( 1 +" < " l) ) : 

k 



— (3 - 2n)(coth(ai)) 2 + 2 h — 



n (siah.(at))« 



A)- 



(m) 



'(sinh(ai))* (1+ " (m>) 

The above solutions satisfy energy conservation equations ((8]) and ((9|) identically. 



(21) 
(22) 
(23) 

(24) 



(25) 



The behavior of EoS for DE (oj de ) in term of cosmic time t is shown in Fig. 2. It is observed that for all three 
closed, flat and open models of the universe, the EoS for DE us de is decreasing function of time, the rapidity of 
their falling down at the early stages depend on the type of the universe, while later on the EoS parameter for 
all three models tend to the same constant value —1 independent to it. We also observe that EoS parameter of 
closed and flat universe are varying in quintessence era (oj de > —1) through out the evolution, while later on 
they tend to the same constant —1 (i.e. cosmological constant) independent to it. We also observe that model of 
open universe started its evolution from quintessence era and crosses the PDL (cj rfe = —1) and ultimately/finally 
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Figure 2: The plot of DE EoS parameter uj de versus t for closed, flat and open universe. Here uj m — O.Of , 
a = 1 .5, n = 1 .5, p$ = 1. 



approaches to — 1 (i.e. cosmological constant). Therefore, we observe that the variation of ui de in our derived 
models is in good agreement with recent observations of SNe la data (Knop et al. [TS] , SNe la data with CMBR 
anisotropy and galaxy clustering statistics (Tegmark et al. [5])- 



The density parameter ft'" 1 ) of perfect fluid and density parameter ft( de ) of DE are given by 

o(m) = P^ = Pon 2 (sinMat))^-^^"""')) 
3H 2 3a 2 (cosh(crf)) 2 



p (de) 

3IP 



3a 2 



(tanh(at)) 5 



(sinh(crf))* (sinh(at))*( 1 +" ( '" ) ) 



Adding ([26]) and (|27| . we get the overall density parameter 



ft = ft (m) + ft (de) = 1 



n tanh(ai) 
a(sinh(at))^ 



(26) 



(27) 



(28) 



From the right hand side of Eq. ([28]) it is clear that in flat universe (k = 0), ft = 1 and in open universe 
(A; = —1), ft < 1 and in closed universe (k — +1), ft > 1. But at late time we see for all flat, open and closed 
universes ft — > 1. This result is also compatible with the observational results. Since our model predicts a flat 
universe for large times and the present-day universe is very close to flat, so the derived model is also compatible 
with the observational results (Bennett et al.[7]; Tegmark et al.[8]). The variation of density parameter with 
cosmic time has been shown in Fig. 3. 



5 Discussions and Concluding Remarks 

In this present work we continue and extend the previous work of Amirhashchi et al. |17j and Saha et al. 
[20] . In summary, we have studied a system of two fluid within the scope of a spatially homogeneous and 
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Figure 3: The plot of ft m , n de , ft versus t for flat universe (k = 0). Here uj m = 0.01, a = 1.5, n = 1.5, p = 1. 

isotropic FRW model. The role of two fluid minimally coupled in the evolution of the dark energy parameter 
has been investigated. The field equations have been solved exactly with suitable physical assumptions. The 
solutions satisfy the energy conservation equation identically. Therefore, exact and physically viable FRW 
model has been obtained. It is to be noted that our method of solving the held equations is different from the 
technic of Kumar [21j . Kumar has solved the held equations by considering the constant DP whereas we have 
considered time-dependent DP. As we have already mentioned in previous section that for a universe which was 
decelerating in past and accelerating at the current epoch, the DP must show signature flipping ( Padmanabhan 
and Roychowdhury [28]; Amendola 29 ; Riess et al. 30 ). So, it is reasonable to consider time dependent DP. 
The main features of the model are as follows: 

• The present DE model has a transition of the universe from the early deceleration phase to the recent 
acceleration phase (see, Figure 1) which is in good agreement with recent observations (Caldwell et al. 
02]). 

• The DP (q) as a function of the red shift parameter z = — 1 + where do is the present value of the scale 
factor i.e. at z = 0, is given by 



q(z) — n — 1 — n 



tanh sinh 



n - 1 - go 
(g + l)(2 + l) 2 " 



(29) 



Here qo is the present value of deceleration parameter i.e. at z — 0. If we set go = —0.73 (Cunha and Lima 
[33] ) for the present Universe (to = 13.7 GYr), we get the following relationship between the constants n 
and a: 



13.7 



■ tanh 



1 - 



0.27 



(30) 



It is self explanatory from the above relation that for the present universe, the model is valid only for 
n > 0.27. Figure 4 depicts the behaviour of q with red shift z for qo = —0.73 and for a representative case, 
we have chosen n = 3 and a — 0.137 satisfying the above Eq. ([30]) . It is clearly observable from the Fig. 4 
that the transition red shift (z t ) for our model with qo = —0.73 is = 0.32. This is in good agreement with 
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Figure 4: The plot of deceleration parameter q vs redshift z. Here a = 0.137, n — 3. 

the cosmological observations in the literature (Cunha and Lima [33]; Cunha [34]; Pandolfi [35] : Lima et 
al. [36]; Li et al. [37]), according to which the transition red shift (zt) of the accelerating expansion is 
given by 0.3 < z t < 0.8. In particular, the kinematic approach to cosmological data analysis provides a 
direct evidence to the present accelerating stage of the universe, which does not depend on the validity of 
general relativity, a well as on the matter-energy content of the universe (Cunha and Lima [33] ) . 

• It is observed that EoS parameter of closed and flat universe are varying in quintessence era (uj de > — 1) 
through out the evolution, while later on they tend to the same constant —1 (i.e. cosmological constant) 
independent to it. 

• It is observed that the open universe started its evolution from quintessence era and crosses the PDL 
(oj de = —1) and finally approaches to —1 (i.e. cosmological constant). Thus, we find that the EoS parameter 
for open universe changes from oj de > — 1 to uj de < — 1, which is consistent with recent observations. 

• The total density parameter (f2) approaches to 1 for sufficiently large time (see, Figure 3) which is repro- 
ducible with current observations. 

• For different choice of n, we can generate a class of DE models in FRW universe. It is observed that such 
DE models are also in good harmony with current observations. For example, if we put n = 1 in the 
present paper, we obtain all results of recent paper of Amirhashchi et al. [17] . 

• Thus, the solutions demonstrated in this paper may be useful for better understanding of the characteristic 
of DE in the evolution of the universe within the framework of FRW space-time. 
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